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A new effective interaction PKA1 with p-tensor couplings for the density-dependent rel- 



j ^ . ativistic Hartrce-Fock (DDRHF) theory is presented. It is obtained by fitting selected cm- 

■ pirical ground state and shell structure properties. It provides satisfactory descriptions of 

nuclear matter and the ground state properties of finite nuclei at the same quantitative level 
as recent DDRHF and RMF models. Significant improvement on the single-particle spectra 
is also found due to the inclusion of p-tensor couplings. As a result, PKAl cures a common 
disease of the existing DDRHF and RMF Lagrangians, namely the artificial shells at 58 



and 92, and recovers the realistic sub-shell closure at 64. Moreover, the proper spin-orbit 
splittings and well-conserved pseudo-spin symmetry are obtained with the new effective in- 
teraction PKAl. Due to the extra binding introduced by the p-tensor correlations, the 



> 

OjS ■ balance between the nuclear attractions and the repulsions is changed and this constitutes 

£C) ' the physical reason for the improvement of the nuclear shell structure. 
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I. INTRODUCTION 

Within the relativistic scheme, mean field theory has achieved great success in the description of 
finite nuclei and nuclear matter during the past years. One of the most outstanding models is the 
relativistic 
theory 

exotic nuclei with a limited number of free parameters, i.e., meson masses and meson-nucleon coupling 

i y, y, q, y, y, 



Hartree approach with the no-sea approximation, namely the relativistic mean field (RMF) 

m 

. The RMF theory provides appropriate quantitative descriptions for both stable and 



constants 



10 
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121 ]. Especially, the RMF model provides a natural mechanism for 
explaining the spin-orbit splittings in nuclear spectra with the covariant formulation of the strong 
scalar and vector fields. This feature becomes even more of central importance with the experimental 



*Electronic address: whlong@pku.org.cn 



observation that nuclei near drip lines undergo a modification of their shell structure, where the 

spin-orbit potential must play an essential role. 

In the framework of the RMF approach, however, there exist two serious defects. One is the missing 

of one-pion exchange process. Namely, with the relativistic Hartree approach, the one-pion exchange 

has zero contribution in mediating nuclear interactions. Because of its small rest mass, one cannot 

simulate the one-pion exchange contributions with zero-range limit. It is expected that the one-pion 

exchange may have minor effects in the spin-saturated system, whereas in the unsaturated systems 

it plays an essential role in determining the isospin dependence of the shell evolutions [0]. Another 

problem is the tensor correlations, e.g., the p-tensor couplings. In the Hartree approximation, the 

contributions from the tensor couplings are practically negligible. In recent non-relativistic 0] 

and relativistic studies UM , it is shown that the tensor forces have distinct effects on the shell evolution 
□ 

of nuclei [16]. In fact these two defects are mainly due to the absence of Fock terms which are dropped 
in RMF. 

During the past decades, there have been several attempts to include the Fock terms in the relativis- 



tic description of nuclear systems 



17 



approaches 

could not provide satisfactory quantitative descriptions for the nuclear structure properties compared 
to the RMF approach. This is mainly due to the numerical complexity induced by the inclusion of Fock 
terms, which strongly increase the difficulties to find appropriate effective Lagrangians for the RHF 
approach. Recently a new RHF method, the density-dependent relativistic Hartree-Fock (DDRHF) 
theory 



221, 



2l| has brought a new insight to this problem. With the effective Lagrangians of Refs. {21 . 
the DDRHF theory can describe the ground state properties of many nuclear systems quantitatively 
on the same level as RMF. In addition, the investigations about the nuclear shell structure evolution 
within the DDRHF theory indicate that the one-pion exchange has a significant effect on the isospin 
dependence of the shell evolution 0]. 

However, some artificial shell structures, e.g., Z = 58 and Z = 92 appear in the calculations of 
RMF 



em in these regions 



24] 



231 ] as well as RHF. These spurious shells lead to the overbinding prob 
and they also affect strongly the isospin dependence of the shell evolutions [13[ . The relative positions 
of l<?7/2 and 2c?5/2 states induce an artificial shell closure at Z = 58. The corresponding states for 
Z = 92 are lhgn and 2/ 7 / 2 - As a common feature, all these states are high-j states. Then, the 
single-particle energies of these states will be strongly affected by the tensor force [0], e.g., p-tensor 
couplings, which was not included in RMF or DDRHF before. In order to solve this artificial shell 
structure problem, we consider the p-tensor correlations in this work. In Section HI we introduce the 
p-tensor couplings in the DDRHF theory, where a new effective interaction PKAl with the p-tensor 
coupling is presented. In Section Hill the detailed investigation of the nuclear structure is performed 
with the newly obtained effective interaction PKAl. The conclusions are drawn in Section IIVI 



3 



II. DDRHF THEORY WITH p-TENSOR COUPLING 



The starting point of the DDRHF theory is the Lagrangian associated with nucleon (tp), isoscalar a 



21 



221 ] . In the isoscalar channels, the 



and (j-mesons, isovector p- and 7r-mesons, and photon (A) fields 
cr-scalar and w- vector couplings provide the main part of the nuclear interactions, i.e., the short-range 
repulsive and mid-, long-range attractive interactions, respectively. One of the distinct differences of 



RHF from RMF is that all the mesons, including t 
the isospin part of nuclear interactions. In Refs. 



re isoscalar ones, have significant contributions to 
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22. 



251 ] . the p- vector and 7r-pseudo- vector 



couplings were introduced in the calculations of DDRHF. The recent investigation about the role of 
one-pion exchange in DDRHF shows that the tensor type force has the strong effects on the nuclear 
structure [l3]. In this study, we introduce the p-tensor correlations into the DDRHF theory in order 
to have a better understanding of the shell structure and cure the artificial shell structure problems 



of DDRHF and RMF 



23] 



A. General Formalism for p-Nucleon Couplings 



In the DDRHF theory, the part of the Lagrangian containing the p-meson fields can be written as 



s? 



i 



i 



[1) 



where R pu = d p p u — d u p p , and m p denotes the rest mass of the isovector-vector p-meson (p v ), and g p 
and f p are the vector and tensor coupling strengths, respectively. 

From the Lagrangian ([T]), one can obtain the equation of motion for the p-meson field as, 



f 

[□ + m 2 p ] p v = gp^yTip+d 11 ^^a U pf^ , 



which leads to the general form of p- meson field, 

p v (x) = / d 4 y 



(2) 



(3) 



where D p is the p-meson propagator. 

With the Legendre transformation {(pi is the field variables of the Lagrangian Jzf), 

dJ? 



J?p 



the Hamiltonian of the p-meson field can be obtained as, 



gptp^Tip—^tpa^ftpdi 



• Pv(x) . 



(4) 



(5) 



In the above expression, we neglect the time component of the four-momentum carried by the mesons, 
which amounts to ignoring the retardation effects. 



With the quantization of the nucleon field ij) [181 ] and the general form of the p- meson field ([3]) , the 



Hamiltonian H p can be written as, 



Hf = -y 



(6) 



x y d ri dr 2 f a ( ri )fp(r 2 ) [r^(l,2>(m p ;l,2)] / 7 (r 2 )/,y(ri) , 
where the Yukawa factor v(m p ; 1, 2) is 



u(m p ; 1,2) 



e ~m P |ri-r 2 



(7) 



47r |ri — r 2 | 

The inclusion of p-tensor correlations leads to three types of interactions in Eq. ([U]): the Vector (V), 
Tensor (T) and Vector- Tensor (VT) couplings. The corresponding vertex matrix T p (l,2) in Eq. (|6j) 
can be expressed as, 



r;(l,2) = ^(1)7^(2)7^(2) 
1 



rj(i,2) 



rY T ( 



1,2) 



p f p (l)a 1/k (l)f p (2)a vl (2)d k (l)d l (2), 

1 



(8b) 



2M L 



^(1)^(1)^(2)7,(2)^(1) + g p (l) lv f p {2)a« v (2)d k {2) \ . 



ku i 



Detailed expressions about the energy functional of p-meson field can be found in Refs. [18l. l21l|. 

In DDRHF, the meson-nucleon coupling constants are treated as a function of baryonic density pb- 
Here we take the same functional form for the density-dependence of the isoscalar mesons (g a and g w ) 



as in Ref. 



2l|: 



9i(Pb) = 9i(po)fi(Q, fovi = (T,uj, 



where 



MO = a 



1 +bi(i + dif 



(9) 



(10) 



1 + + 

with ^ = Pb/po, and po denotes the saturation density of nuclear matter. For the isovector mesons as 
well as the newly introduced p-tensor coupling f p , the exponential density-dependence is adopted as, 



9i = 9i(0)e- a ' * 



(11) 



In the above expression, gi(0) corresponds to the free coupling constants g p , f p , and and cij are the 
corresponding parameters a p , ax and a n , respectively. 

For the open shell nuclei, the pairing correlations are treated by the BCS method and the pairing 



matrix elements are calculated with a zero-range, density-dependent interaction 



26] 



V(r 1 ,r 2 ) = V 5(ri - r 2 ) 



Pb{r) 



(12) 



5 

where V = -900 MeV-fm 3 . The active pairing space is limited to the single-particle states below the 
single-particle energy +15 MeV. 

In this work, the corrections from the center-of-mass motion are treated in the same way as in 



Refs. 
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221 ] . For the numerical calculations, a box boundary condition at 20fm is introduced for the 
unbound states as well as the bound ones and we check that the overall results are not affected by 
the choice of the box size. For the radial step, one may choose smaller one about 0.05 fm in the light 
nuclei whereas 0.1 fm is precise enough for the heavy nuclei. 



B. New Effective Interaction 



In the previous DDRHF parametrizations (PKOl, PK02 and PK03) [lj, |2l|, |22j , we selected 12 
nuclei as the reference ones, i.e., 16 0, 40 Ca, 48 Ca, 56 Ni, 68 Ni, 90 Zr, 116 Sn, 132 Sn, 182 Pb, 194 Pb, 208 Pb 
and 214 Pb. In the present case for PKAl, we aim to cure the artificial shell structure problem at 
Z = 58 and 92. To solve this problem and recover the sub-shell closure at Z = 64, two more nuclei 
14 gCe82 and 14 |Gds2 are added as the reference nuclei and 56 Ni is replaced by its neighboring one 



11 



22|. Besides the bulk properties 



58 Ni. The parameter fitting procedure is similar to that in Refs. 
(poj K and J) of nuclear matter and the binding energies of the reference nuclei, we include the 
spin-orbit splittings of neutron and proton Ip states of 16 0, and the shell gaps at Z = 58 ( 14 gCe82) 



22], we obtain 



and Z = 64 ( 14 f Gdg2) as the new criteria. By minimizing the x 2 error as in Ref. 
the new effective interaction PKAl with the p-tensor coupling for the DDRHF theory (see Table Q]). 
In this parametrization, we have 12 free parameters, 6 in the isoscalar channels as well as 6 in the 
isovector channels. We slightly change the coupling strength for 7r-meson (/,,- and a n ) in PKAl from 



the effective interaction PKOl 



211 ] . which is the starting Lagrangian in the fitting process. 



TABLE I: Parameters of new effective interaction PKAl. The quantity k is k = f p /g p , the ratio of the vector 
and tensor couplings of p-meson. The masses (in MeV) of nucleon, uj-, p- and 7r-mesons are taken as M = 938.9, 
to w — 783.0, nip = 769.0 and = 138.0, respectively. 





488.227904 


a a 


1.103589 




1.126166 




8.372672 


b a 


16.490109 




0.108010 


9u 


11.270457 


Ccr 


18.278714 




0.141251 


Po 


0.159996 


d a 


0.135041 




1.536183 


9 P 


3.649857 


K 


3.199491 




1.030722 


a p 


0.544017 


ax 


0.820583 




1.200000 



In Fig. dj the density-dependent couplings g a , g u (left panels) and g p , / p , and f n (right panels) 
of the effective interaction PKAl are shown as functions of the baryonic density p^, in comparison 
with PKOl (DDRHF) 21] and DD-ME2 (RMF) Q- As seen from Fig. [TJ the density-dependence 



6 

of isoscalar couplings, especially for g a , is weak for PKAl compared to PK01 and DD-ME2. For 
the isovector channels, g p , f„ and f v show strong density-dependence in PKAl whereas the density- 
dependence of g p of PKOl is weak. For the one-pion exchange, PKAl and PKOl have nearly the 
same coupling strength for f w . PKAl and PKOl have smaller coupling strength in magnitude for 
both the isoscalar and isovector channels compared to DD-ME2. This is mainly due to the effects 
of the Fock terms. One can also find that the inclusion of />tensor coupling leads to smaller g a and 
larger g w in PKAl than those in PKOl. This indicates that the p-tensor correlations contribute to 
make the nuclear interactions attractive. 
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FIG. 1: (color online) The density-dependent mcson-nucleon couplings in the isoscalar (left panel: g a and 
gui) and isovector (right panel: g p , f p and f^) channels as functions of density for the new DDRHF effective 
interaction PKAl, in comparison with PKOl in DDRHF and DD-ME2 in RMF. The shadowed area denotes 
the empirical saturation density region. 



We calculate the bulk properties of nuclear matter with the effective interaction PKAl as shown in 
Table [Ml where the results calculated by DDRHF with PKOl and RMF with DD-ME2 are also listed 
for comparison. Compared to PKOl and DD-ME2, PKAl gives a larger saturation density po, which 
is close to the common value in non-relativistic HF calculations. Although PKAl gives a smaller 
compression modulus K and a larger symmetry energy J, the values are still acceptable. Among 
these three effective interactions, PKOl has the largest effective masses (the relativistic one and 
the non-relativistic one M^ R ) [21( whereas DD-ME2 gives smallest ones. Comparing the values of 
and M%, one can find a significant difference between DDRHF and RMF, the DDRHF models 
giving a larger difference between these two masses. Actually this difference is of special importance 



7 

in describing nuclear structure since the effective mass Mj^ R is related to the level density whereas the 
scalar mass M$ is related to the spin-orbit splitting. In the following, one may find the corresponding 
effects due to this difference. 



TABLE II: The saturation density po (fm -3 ), the binding energy Eb (MeV), the compression modulus K 
(MeV), the symmetry energy J (MeV), the effective masses and M^ R 2l|, and the scalar mass Mg for the 
symmetric nuclear matter. The results are calculated with PKAl and PKOl in DDRHF, and with DD-ME2 in 
RMF. 





Po 


Eb 


K 


J 


Mr 


^NR 


Ms 


PKAl 


0.160 


-15.83 


229.96 


36.02 


0.663 


0.681 


0.547 


PKOl 


0.152 


-16.00 


250.24 


34.37 


0.727 


0.746 


0.590 


DD-ME2 


0.152 


-16.14 


250.97 


32.30 


0.635 


0.652 


0.572 



As we have mentioned before, 14 reference nuclei are adopted in the parametrization of PKAl. 



The binding energies and charge radii of these nuclei calculated with PKAl are shown in Tab 
comparison with the calculations of PKOl 



28 



eO 



3 



in 



21J, DD-ME2 [27|] , and the experimental data 
Other reference nuclei than the present 14 ones used in the parametrizations PKOl and DD-ME2 are 
listed in the lower panel of Table lllll In this table, the root mean square deviations (rmsd) A and 
relative rmsd 5 are defined as 



1 N 

N ^ 

i=l 

N 



Exp. 



N 



y 



Cal. 



V, 



CaL /y E*p. 



(13a) 
(13b) 



From the values of A and 5 in Table lllll we can see that PKAl provides an appropriate description 
for both binding energies and charge radii of these reference nuclei. For the binding energies, PKAl 
gives the best agreement with the data. For the charge radii, PKAl gives quantitatively comparable 
description to PKOl and DD-ME2. Since the reference nuclei cover from light ( 16 0) to heavy ( 214 Pb, 
210 Po) ones, one may expect that DDRHF with the tensor interaction PKAl can provide a proper 
description of the nuclei in the whole nuclear chart. 

The charge densities calculated with PKAl and PKOl for the nuclei 16 0, 40 Ca, 48 Ca, 90 Zr and 
208 Pb are presented in Fig. [2l The experimental data 29[] are also shown for comparison. One can 
see that DDRHF with PKAl provides a fairly good agreement with the data in heavy systems, e.g., 
208 Pb and 90 Zr. For 40 Ca and 48 Ca PKAl also shows comparable quality to PKOl whereas it presents 
less good agreement for 16 0. In Table [TTfl one can also find better agreement in the heavy nuclei than 
in the light ones for the charge radii calculated with PKAl. 



8 



TABLE III: Binding energies and charge radii for the reference nuclei. The results arc calculated by DDRHF 
with the new effective interaction PKA1 and with PKOl [21j, and by RMF with DD-ME2 [27|. The upper 



panel gives the reference nuclei of PKAl and the other reference nuclei of PKOl and DD-ME2 are listed in the 
lower panel. Experimental data arc taken from Refs. 28 
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-126.9913 
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-342.0520 
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3.4423 
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48 Ca 


-415.9904 


-416.3697 
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-505.9705 
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182 Pb 
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194 Pb 
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III. NUCLEAR STRUCTURE PROPERTIES WITH TENSOR CORRELATIONS 

A. Single-particle Spectra 

In Section [TT] we showed that the effective interaction PKAl can describe quantitatively well the 
bulk properties of nuclear matter and the ground state properties of finite nuclei. In this section, we 
study the shell structure of the Hartree-Fock single-particle energies for several reference nuclei by 
using PKAl. For comparison, we also show the results of the reference nuclei calculated with three 



other interactions: PKOl (DDRHF without p-tensor couplings), PK1 11] (RMF with non-linear 
self-couplings of mesons), DD-ME2 (RMF with density-dependent meson- nucleoli couplings). 

In Fig. [3]and Fig. [4]are shown the neutron (left panel) and proton (right panel) single-particle levels 
in 140 Ce and 146 Gd, respectively. For these two N = 82 isotones, their proton numbers correspond 
to the artificial shell closure Z = 58 occurring in RMF and the sub-shell closure Z = 64 observed in 
experiments, respectively. In the results calculated with PKOl, PK1, and DD-ME2, one can find large 
shell gaps at Z = 58 or N = 58 in both nuclei. These gaps are even comparable to the well established 
shell gaps Z = 50 and N = 82. In contrast, PKAl brings much smaller shell gaps at iV = 58 and 
Z = 58 in the single-particle levels of 140 Ce In 146 Gd, the artificial shell structures Z = 58 and N = 58 
completely disappear and the sub-shell closures Z = 64 and iV = 64 are well reproduced by PKAl. 
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FIG. 3: Single-particle energies of 140 Ce. The results are calculated by DDRHF with PKAl and PKOl, and 
RMF with PK1 and DD-ME2. 



Besides 140 Ce and 146 Gd, we also calculated another N = 82 isotone, the doubly magic nucleus 
132 Sn. In Fig. [5] are shown the neutron (left panel) and proton (right panel) single-particle levels 
calculated with PKAl, PKOl, PK1, and DD-ME2. The experimental data from Ref. [3l| are also 
shown for comparison. In this figure, one can also find the spurious shell structures Z or iV = 58 



appearing in the results of PKOl, PK1, and DD-ME2. Compared to these three results, PKAl shows 
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FIG. 4: Single-particle energies of 146 Gd. The results are calculated by DDRHF with PKA1 and PKOI, and 
RMF with PK1 and DD-ME2. 



distinct improvement on this problem. Namely, the spurious shell gaps do not appear any more and 
fairly good agreement with the data is obtained with PKAl. 
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FIG. 5: Single-particle energies of 132 Sn. The results are calculated by DDRHF with PKAl and PKOI, and 
RMF with PK1 and DD-ME2. Experimental data arc taken from Rcf. 0. 



It is well known that the two-nucleon separation energy is an important criteria to identify the 
nuclear shell structure. In Fig. O the two-proton separation energy S*2 P calculated with PKAl, 
PKOI and DD-ME2 are shown as a function of proton number Z for the N = 82 isotones, and the 



experimental data 28| are also shown for comparison. A sudden change of the slope is found in 
the results of PKOI and DD-ME2 at Z = 58. This is a clear sign of the existence of spurious shell 
structure at Z = 58 in the calculations of the existing DDRHF and RMF Lagrangians. In contrast, the 
two-proton separation energy calculated by PKAl decreases smoothly as a function of Z and agrees 
well with the experimental data. Fig. [5] shows that the spurious shell structure Z = 58 is successfully 
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eliminated by the new DDRHF effective interaction PKAl. 

T 




FIG. 6: (color online) Two-proton separation energies of the N = 82 isotones. The results are calculated by 



DDRHF with PKAl and PKOl, and RMF with DD-ME2. The experimental data are taken from Ref. 



28] 



Besides N (or Z) = 58, the nucleon number 92 is another spurious shell structure in the calculations 



of RMF 



231 ] . In Fig. [7] it is clearly shown in both neutron and proton single-particle spectra that 



iV (or Z) = 92 becomes a spurious shell structure in the results of PKOl, PK1, and DD-ME2. With 
the inclusion of the p-tensor correlations, the artificial shell closures at A = 92 and Z = 92 disappear 
in both neutron and proton spectra in the results of PKAl. In addition, PKAl also shows another 
improvement in the order of single-particle levels, e.g., for the neutron states 2<7g/ 2 and Hxi/2 111 208 Pb. 
Namely, DDRHF with PKAl provides the same ordering with the data while the other three cases 
fail to predict the correct ordering. 
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FIG. 7: Single-particle energies of 208 Pb. The results are calculated by DDRHF with PKAl and PKOl, and 
RMF with PK1 and DD-ME2. The experimental data are taken from Ref. 



13 



B. Spin-orbit Splittings and p- Tensor Correlations 

From the discussions on the single-particle spectra of several reference nuclei, one can find signifi- 
cant improvements with the inclusion of p-tensor correlations in DDRHF calculations. The spin-orbit 
splitting is also very essential to test the validity of the model. In Table HVi we give the spin-orbit 
splittings in several magic nuclei obtained by using PKAl, PKOl, and DD-ME2. The experimental 
data are also tabulated for comparison. As seen from this table, PKAl provides comparable quanti- 
tative results of the spin-orbit splittings to PKOl and DD-ME2. Among the results calculated with 
PKAl, there are some systematic over estimations of the spin-orbit splittings of some states, e.g., z^lf 
in 48 Ca and 56 Ni, z^lg in 90 Zr, 7iTg in 132 Sn, and z^li, 7rlh in 208 Pb, which account for the corresponding 
shell closures. For these spin partner states, one of the partners is not occupied and one may expect 
sizable corrections from particle-vibration coupling. In the present calculations, the particle- vibration 
coupling is not included yet, and its effects generally tend to shift the occupied and unoccupied states 
to the Fermi surface. Therefore, the systematic over estimation of the spin-orbit splittings may leave 
some space for this effect. 



TABLE IV: The spin-orbit splittings (in MeV) of neutron (y) and proton (tt) states for the magic nuclei 
calculated with PKAl, PKOl, and DD-ME2. The experimental data are taken from Rcf. 31|. 



Nucleus 


State 


Exp. 


PKAl 


PKOl 


DD-ME2 


Nucleus 


State 


Exp. 


PKAl 


PKOl 


DD-ME2 


16Q 


z^lp 


6.18 


6.055 


6.426 


6.545 


56 Ni 


vli 


6.82 


10.027 


7.363 


8.361 




7rlp 


6.32 


5.973 


6.356 


6.472 




v2p 


1.11 


0.899 


0.977 


1.383 


40 Ca 


z/ld 


6.75 


7.386 


6.742 


6.760 


90 Zr 


v2p 


0.37 


1.702 


1.598 


1.686 




u2p 


2.00 


2.527 


1.846 


1.694 




v2d 


2.42 


2.457 


2.049 


2.097 




TTld 


5.94 


7.215 


6.629 


6.696 




vlg 


7.07 


8.608 


7.155 


7.609 


48 Ca 


vld 


5.30 


6.817 


5.414 


6.172 




tt2 P 


1.51 


1.489 


1.430 


1.619 




ulf 


8.01 


8.512 


7.345 


7.737 


208p b 


^2f 


2.14 


2.342 


2.009 


2.317 




v2p 


1.67 


1.647 


1.347 


1.462 




is2g 


2.38 


2.482 


2.103 


2.322 




TTld 


5.01 


6.833 


5.590 


6.406 




vli 


5.81 


7.936 


6.143 


6.970 




?r2p 


2.14 


1.634 


1.322 


1.539 




i/3p 


0.90 


0.925 


0.782 


0.889 


132 Sn 


v2d 


1.66 


1.866 


1.645 


1.969 




i/3d 


0.89 


0.875 


0.710 


0.726 




Trig 


5.33 


6.909 


5.291 


6.438 




7r2d 


1.53 


1.506 


1.404 


1.769 




Tr2d 


1.75 


1.569 


1.463 


1.912 




Trlh 


5.03 


6.380 


4.901 


5.955 



From the Dirac equation, one can express the single-particle energy for a state a as 
E a = Ek,a + E(T,a + E UA + E P:a + E Wja + Ea,cl + E^ a , 



(14) 



where E^ ^ denotes the kinetic contribution, and Ei A (i = a,uj, p,ir, A) represent the contributions 
from the mesons and. photon coupling channels including the direct and exchange parts, and E^ a 
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accounts for the rearrangement terms. From Eq. (|14p . one can also obtain the contributions to the 
spin-orbit splittings from different channels. 
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FIG. 8: Contributions from p- and 7r-mesons to the spin-orbit splittings for the neutron orbits in 146 Gd calculated 
by DDRHF with PKA1 and PKOI. The left panel shows the total values of the spin-orbit splittings and the 
contributions from the p- and 7r-mcsons, and the right panel shows the ratios of the p- and 7r-mcsons contributions 
to the spin-orbit splittings. 
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FIG. 9: Similar as Fig. El but for 208 Pb. 



In order to understand the improvement on the shell structure brought by the p-tensor couplings, 
we compare the contributions of the p- and 7r-mesons to the spin-orbit splittings between two DDRHF 
effective interactions PKAl and PKOI. In Fig. [B] and Fig. [5J are shown the contributions to the 
spin-orbit splittings from p- and ir- mesons, respectively for the 146 Gd and 208 Pb neutron orbits. From 
these two figures, one can find that the tt pseudo- vector (PV) coupling gives almost same contributions 
to the spin-orbit splittings in PKAl (filled up-triangles) and PKOI (open up-triangles) . This can 
be well explained by the equivalent 7r-coupling strength in these two effective interactions. For the 
contributions from /9-meson, PKAl and PKOI have distinct difference. Due to the tensor interactions, 
the p- meson coupling in PKAl (filled circles) shows substantially larger effects than that of PKOI 
(open circles). It is seen that the average contributions from the p-meson couplings in PKAl are 
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about 20% of the total ones. For the states near the Fermi surface or some high-j orbits, the p-meson 
contribution grows up to 30% and even higher. On the other hand, the p-meson coupling of PKOl 
gives about 10% in most cases. The 7r-meson gives opposite contributions to the p-meson for the spin- 
orbit splittings except the lp states and the magnitude is about (10~ 20)% of the spin-orbit splitting. 
Because of the tensor effects in 7T-PV couplings, some systematic enhancements are also observed in 
high-j states like lg and lh states in 146 Gd, and lh and 1% in 208 Pb. As was noticed in the comparison 
between PKAl and PKOl, the p-tensor couplings have significant effects on the spin-orbit splittings, 
especially for high-j orbits, and affect much the shell structures. This is the main reason why the 
improvement of the shell structure is obtained with the />tensor correlations. 



C. Spurious Shell closures and Pseudo-spin symmetry 

As we mentioned before, the spurious shell closures at 58 and 92 are related with the pairs of high-j 
states {2^5/2,157/2} and {2/7/2,1/19/2}; respectively. These pairs are the pseudo-spin partners, If 



and lg states, respectively. The spurious she 
pseudo-spin symmetry (PSS) 
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32, 



33, 



34, 



closure problem is then related to the conservation of 



36. 



371 ] . i.e., the existing artificial shell structures in 



RMF break largely PSS. As seen from the results of PKAl (see Fig. [4] and Fig. [7]), PSS is successfully 
recovered for the If states in 146 Gd and lg states in 208 Pb. In order to understand the improvement 
due to the p-tensor correlations, we studied the contributions from different terms in Eq. (|14p to the 
pseudo-spin orbital splittings. In Table [V] and Table [\7TI are shown the results calculated by DDRHF 
with PKAl (upper panels) and PKOl (lower panels) respectively for 146 Gd and 208 Pb neutron orbits. 
One can find in these results that PKAl conserves pseudo-spin symmetry better than PKOl for the 
states near the Fermi levels, e.g., vlf and u2p states in 146 Gd, vlg and vld states in 208 Pb. 

For the pseudo-spin orbital splittings, PKAl and PKOl provide similar contributions in magnitude 
to the kinetic part (AE^), the rearrangement term (AEr), and the 7r-coupling (AE n ) except for a few 
cases. For the contributions from a-, w-mesons (AE a+UJ ) and p-meson (AE p ) couplings, there exist a 
distinct difference between PKAl and PKOl, especially for the states near the Fermi surfaces. From 
these two tables one can see that the p-meson couplings in PKAl give larger contributions to the 
pseudo-spin orbital splittings than those in PKOl and the p-tensor couplings increase the splittings. 
For the states near the Fermi surface, PKAl provides negative values of AE a+UJ , which cancel largely 
with AEk and AE p . In the PKOl results, the AE a+UJ is always positive and only the rearrangement 
term AEr partially cancels the contributions from the other channels. 

To understand the differences between the results of PKAl and PKOl, we study the contributions 
from different terms to the average binding energy E of the spin partner states ji and j2, i.e., E = 
[Ei(2ji + 1) + E 2 (2j 2 + 1)] / (2ji + 2j 2 + 2). In the left panel of Fig. [101 the average binding energies 
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TABLE V: The contributions (in MeV) from different terms in Eq. (fT4|) to the pseudo-spin orbital 
splittings AE for neutron (v) orbits in 146 Gd , calculated by DDRHF with PKA1 (upper panel) and 
PKOl (lower panel). The average binding energy E for the pseudo-spin partner states j\ and ji is 
[£7 1 (2j 1 + 1) + E 2 (2 ]2 + 1)] / (2 :h + 2j 2 + 2). 



146 Gd 


State 


E AE 


AE P AE n AE R AE k AE a+U} 


AE a AE^ 


PKAf 


vlp 
vld 

vlf 

v2p 


-39.44 2.935 
-27.68 2.097 
-16.01 0.489 
-13.35 0.422 


0.277 0.039 0.080 1.900 0.639 
0.476 0.225 -0.439 2.592 -0.756 
1.107 0.525 -1.187 2.546 -2.502 
0.602 0.344 -0.456 0.854 -0.920 


11.391 -10.751 
14.847 -15.603 
18.145 -20.647 
7.499 -8.419 


PKOl 


vlp 
vld 

fl/ 

v2p 


-37.63 2.985 
-27.02 3.013 
-16.13 2.224 
-12.68 0.851 


0.429 0.041 -0.370 1.521 1.364 
0.493 0.180 -0.814 2.131 1.023 
0.662 0.404 -1.360 1.945 0.573 
0.387 0.284 -0.524 0.188 0.516 


8.502 -7.138 
12.372 -11.349 
18.306 -17.733 

9.315 -8.799 



TABLE VI: Same as Table for 208 Pb neutron orbits. 



3 Pb 



State 



E 



AE 



AE n AE„ AE 



R 



AE k AE, 



<7-\-LO 



AE n 



AE„. 



PKA1 



vlp 
vld 

"If 
vlg 

v2p 

v2d 



-42.74 
-32.46 
-22.06 
-11.87 
-18.76 
-9.17 



3.292 
3.560 
2.502 
0.584 
0.259 
0.092 



0.342 
0.329 
0.608 
1.181 
0.265 
0.777 



0.083 
0.218 
0.455 
0.757 
0.233 
0.459 



-0.006 
-0.111 
-0.563 
-1.251 
-0.399 
-0.695 



1.655 
2.393 
2.738 
2.431 
1.177 
0.773 



1.218 
0.730 
-0.735 
-2.535 
-1.017 
-1.222 



7.555 
14.699 
20.507 
23.240 
3.704 
5.289 



-6.337 
-13.969 
-21.242 
-25.776 
-4.721 
-6.510 



PKOl 



vlp 
vld 

vlg 

v2p 
v2d 



-40.04 
-31.04 
-21.66 
-12.18 
-17.95 
-8.67 



2.851 
3.644 
3.397 
2.314 
0.678 
0.547 



0.238 
0.421 
0.549 
0.619 
0.117 
0.183 



0.100 
0.207 
0.380 
0.604 
0.181 
0.371 



-0.349 
-0.712 
-1.125 
-1.530 
-0.348 
-0.566 



1.370 
2.009 
2.212 
1.731 
0.624 
0.124 



1.492 
1.720 
1.381 
0.890 
0.105 
0.434 



4.663 
10.406 
16.130 
21.939 
4.354 
8.938 



-3.172 
-8.686 
-14.749 
-21.048 
-4.250 
-8.503 



E and the sum of the kinetic part (E^), the a- and w-coup lings (E a and EJ) are shown as a function 
of angular momentum I for the neutron orbits of 208 Pb. The right panel shows the rearrangement 
term (En), and the sum of the p- and 7r-coupling terms {E p and E n ). The results are calculated with 
PKAl (filled symbols) and PKOl (open symbols). In the left panel, it is found that PKA1 provides 
stronger Z-dependence than PKOl for the average binding energy E. In the results of PKOl, the 
main contribution to E is given by the sum Ek + E a + E^, whereas in the results of PKAl, Er and 
E p + En also provide significant contributions to E. This is due to the fact that Er and E p + E n 
cancel each other in the results of PKOl as shown in the right panel. On the other hand, PKAl gives 
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weaker rearrangement term En, and much stronger p-couplings than PK01. Thus, the inclusion of 
/>tensor couplings give significant contributions to the nuclear attraction, which strongly affects on 
the coupling strength in other channels, e.g., PKAl has a stronger ^-coupling than PKOl as shown 
in Fig. d 
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FIG. 10: The average binding energy E = [£'i(2ji + 1) + E 2 (2j 2 + 1)] / {2j 1 + 2j 2 + 2) of the 208 Pb neutron 
orbits as a function of the angular momentum I calculated by DDRHF with PKAl (filled symbols) and PKOl 
(open symbols). Left panel gives the sum contributions from the kinetic part (E^), the a- and w-couplings {E a 
and E u ) to the average binding energy and the total ones, and right panel shows the contributions from the 
rearrangement term (En), the p- 1 and 7r-couplings (E p and E n ). 



In Fig. 1111 are shown the values of E a and E u calculated with PKAl (filled symbols) and PKOl 
(open symbols) as a function of angular momentum /. One can see that PKAl leads to a stronger 
/-dependence of E w than PKOl, and to a similar (slightly stronger) /-dependence of E a for the states 
near the Fermi surface. It should be noticed that the pseudo-spin partner states have different angular 
momenta /. The stronger /-dependence of E u given by PKAl leads to larger negative contributions 
to the pseudo-spin orbital splittings as shown in the last column of Table [V] and Table I VII These 
results finally induce negative values for AE a+w so that PSS can be well conserved in the results of 
PKAl. As a relativistic symmetry, the conservation of the PSS is mainly determined by the balance 
of the nuclear attractions and repulsions 



3a ]. which is also well demonstrated by Table [V] and Table 

IVI1 Compared to the PKOl results, this balance is much changed by PKAl due to the extra binding 
induced by the p-tensor couplings, which indicates the physical reason for the improvement of the 
nuclear shell structure. 
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FIG. 11: Contributions from the a- and w-couplings (E a and E u ) to the average binding energy E = 
\E\(2ji + 1) + i?2(2j2 + 1)] / (2ji + 2j-2 + 2) of the 208 Pb neutron orbits as a function of the angular momentum 
I calculated with PKAl (filled symbols) and PKOl (open symbols). 

IV. CONCLUSIONS 

In this work, we have introduced the />tensor correlations in the density-dependent relativistic 
Hartree-Fock (DDRHF) theory. By fitting the empirical properties of ground state and the shell 
structure, we propose a new DDRHF effective interaction with p-tensor couplings, PKAl. With the 
newly obtained effective interaction PKAl, DDRHF provides satisfactory descriptions of the bulk 
properties of nuclear matter and the ground state properties of finite nuclei, at the same quantitative 
level as the established DDRHF and RMF models. 

Moreover, the inclusion of p-tensor correlations brings a significant improvement on the descriptions 
of nuclear shell structures compared to the existing DDRHF and RMF Lagrangians. Particularly, 
we have studied the sing le-particle spectra of nuclei 140 Ce, 146 Gd, 132 Sn and 208 Pb with PKAl and 
compared to previous DDRHF and RMF approaches. It has been found that the previous DDRHF and 
RMF calculations give the spurious shell closures 58 and 92, whereas the realistic sub-shell closure 64 
cannot be well reproduced. The effective interaction PKAl cures these common diseases, eliminating 
the spurious shell structure and recovering the sub-shell closure 64. In addition, the inclusion of tensor 
correlations improves the descriptions of the ordering of the single-particle levels, e.g., the neutron 
states 2g 9 / 2 and lin/2 in 208 Pb, which are important states for nuclear structure problems. 

The spin-orbit splittings and the pseudo-spin orbital splittings of the magic nuclei are also studied 
by using PKAl, and the PKOl version which has no p-tensor coupling. It is shown that the p-tensor 
correlations have substantial effects on enlarging both splittings, especially for the high-j states. Even 
though, PKAl still provides an appropriate quantitative agreement with the experimental data on 
the spin-orbit splittings for the magic nuclei at the same level as the modern DDRHF and RMF 
Lagrangians. It is shown that the artificial shell structure problem is intimately related to the con- 
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servation of the pseudo-spin symmetry, which is determined by the balance of the nuclear attractions 
from ex-meson and p-tensor couplings and the repulsion from w-meson coupling. It is found that a 
better conserved pseudo-spin symmetry is obtained with PKAl, in which the />tensor correlations 
contribute significantly to the nuclear attraction. Due to the extra binding introduced by the p-tensor 
correlations, the balance of attraction and repulsion is changed by the parametrization PKAl, and 
this constitutes the physical reason for the improvement of the nuclear shell structure. 
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